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Abstract. In this paper, we construct the bihnear identities for the wave func- 
tions of an extended Kadomtsev-Petviashvih (KP) hierarchy, which is the KP 
hierarchy with particular extended flows (2008, Phys. Lett. A, 372: 3819). By 
introducing an auxiliary parameter (denoted by z), whose flow corresponds to 
the so-called squared eigenfunction symmetry of KP hierarchy, we find the tau- 
function for this extended KP hierarchy. It is shown that the bilinear identities 
will generate all the Hirota's bilinear equations for the zero-curvature forms 
of the extended KP hierarchy, which includes two types of KP equation with 
self-consistent sources (KPSCS). It seems that the Hirota's bilinear equations 
obtained in this paper for KPSCS are in a simpler form by comparing with 
the results by Hu and Wang (2007, Inverse Problems, 23: 1433). 
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1 Introduction 

Sato theory has fundamental importance in the study of integrable systems 
(see [6] and references therein). It reveals the infinite dimensional Grassman- 
nian structure of space of r-functions, where the r-functions are solutions for 
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the Hirota's bilinear form of Kadomtsev-Petviashvili (KP) hierarchy. The key 
point to this important discovery is a bihnear residue identity for wave func- 
tions called bilinear identity. Bilinear identity plays an important role in the 
proof of existence for r-function. It also serves as the generating function of 
the Hirota's bilinear equations for KP hierarchy |31 |28l [38]. In this paper, we 
will construct the bilinear identity for the extended KP hierarchy [23] . 

The study of integrable generalization is an important subject in the study 
of integrable systems. Several approaches for the generalizations have been 
developed, e.g., constructing new flows to extend original systems. There are 
many ways to introduce new flows to make a new compatible integrable system. 
In [H], the KP hierarchy is extended by properly combining additional flows. 
In [16] , extension of KP hierarchy is formulated by introducing fractional order 
pseudo-differential operators. In [HI [9] , Dimakis and MuUer-Hoissen extended 
the Moyal-deformed hierarchies by including additional evolution equations 
with respect to the deformation parameters. In [3], Carlet, Dubrovin and 
Zhang defined a logarithm of the difference Lax operator and got the extended 
(2+l)D Toda lattice hierarchy. Later, the Hirota's bilinear formalism and 
the relations of extended (2+l)D Toda lattice hierarchy and extended ID 
Toda lattice hierarchy have been studied [3ll[39]. In [19], Li, et al, studied 
the r-functions and bilinear identities for the extended bi-graded Toda lattice 
hierarchy. 

Recently, we proposed a different and general approach to extend (2 + 1)- 
dimensional integrable hierarchies by using the symmetries generating func- 
tions (or the squared eigenfunction symmetries) [25]. This kind of extended 
KP hierarchy can be thought as the generalizations of the KP hierarchy by 
squared eigenfunction symmetries (or ghost flows) [H [351 [361 [37] . The reason 
for constructing such kind of extended KP hierarchy is that it includes two 
types of KP equation with self-consistent sources (KPSCS-I and KPSCS-II), 
which has important applications in hydrodynamics, plasma physics and solid 
state physics (see, e.g., [lOl [III [11 gQl [HI [301 [3ll [321 [33l [131 [11]). The ex- 
tended KP hierarchy also includes the well-known /c-constrained KP hierarchy 
[SI [171 [18] and several (l+l)-dimensional solition equations with sources as re- 
ductions. It has been shown that many (2+l)-dimensional integrable systems 
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can be extended in this way [151 [Ml [291 [12] • We also proposed a gen- 
eralized dressing approach to construct Wronskian type solutions (including 
soliton solutions) for the extended hierarchies [151 [Ml [261 [23 W2\ . 

KPSCS-I and KPSCS-II can also be written in Hirota bilinear forms. In 
[T3l [^. Hu, Wang, et al, proposed a source generation procedure to construct 
the Hirota bilinear form for KPSCS (I and II) on the basis of the bilinear form 
for the original KP equation. 

However, it turns out that it is not an easy task to find Hirota bilinear 
form for KPSCS (I and II). More generally, it remains unsolved to give the 
Hirota bilinear form for every zero-curvature equation in the extended KP 
hierarchy. The most natural way to solve this problem is to consider the 
bilinear identities of the extended KP hierarchy. Because of the importance of 
the bilinear identities in Sato theory, the existence of bilinear identities could 
be an important open question for the extended KP hierarchy. 

In this paper, we will construct the bilinear identities for the extended KP 
hierarchy. Moreover, it will serve as a generating function of all the Hirota 
bilinear forms for the zero-curvature equations in the extended KP hierarchy. 
In particular, it generates Hirota bilinear forms for KPSCS (both type I and II). 
The squared eigenfunction symmetry plays a crucial role in our construction. 
It seems that the Hirota's bilinear equations obtained in this paper for KPSCS 
are in a simpler form by comparing with the results by Hu and Wang |14j . 

This paper is organized as follows. In Section [2l the bilinear identities 
for KP hierarchy with a squared eigenfunction symmetry are constructed. In 
Section [3l by making the squared eigenfunction symmetry as an auxiliary 
flow, we construct the bilinear identities for the extended KP hierarchy. In 
addition, we prove any wave function satisfying bilinear identities will be a 
wave function for the extended KP hierarchy. In Section [H we find the r- 
function for the extended KP hierarchy. We also find the generation functions 
for Hirota bilinear form for the extended KP hierarchy. Subsequently, several 
examples including the Hirota bilinear form for the KPSCS (I and II) and a 
higher order system in the hierarchy are given. In Section [5l we show how to 
go back from Hirota bilinear forms to nonlinear partial differential equations 
(PDEs) for KPSCS, which verifies the correctness of our construction. In the 
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last section, we will give conclusion and remarks, and discuss some problems 
for further exploration. 

2 Bilinear Identities for KP hierarchy with 
the Squared Eigenfunction Symmetry 

Consider the system given in [37], defined by assuming that L = 0+"^^^ Uid~^ 
satisfies both the ordinary KP flows 

dt^L=[Ll,L], (la) 

and a new S^-flow given by 

d,L=[qd~\L]. (lb) 

Here the Lax operator L is a pseudo-differential operator (see [7] for a detailed 
introduction). The "+" sign in subscript part of indicates the projection 
to the non-negative part of with respect to the power of d. There are a 
series of time variables tn, and we can set x = ti according to ( ITal) . 

A sufficient condition for the commutativity of dz- and 9t„-fiow is obtained 
by imposing q and r to be an eigenfunction and an adjoint eigenfunction of 
KP hierarchy, respectively. Namely, we should assume 

dt,^q = Lliq), (Ic) 
a,„r = -{LlYir), (Id) 

for all G N. The S^-fiow describes a symmetry for KP hierarchy, which is 
called the squared eigenfunction symmetry [37] or a ghost flow [1] [2] . 

It is well-known that dressing operator W = J2i>o'^i^~^ (""^o = 1) plays a 
very important role in KP theory. Here Wi {i > 0) are functions of all t„ and 
z. Suppose the t„-evolution equations for W are given by 

dt^W = - {Wd''W-^)_ W, (n G N). (2) 
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Then Lax operator L is related to the dressing operator W hy L = WdW ^. 
We define the wave and adjoint wave function as 

M;(z,t,A) = iye«(*'^\ (3a) 
M;*(z,t,A) = (iy*)"^e-«(*'^) (3b) 

where ^(t, A) = X^iM^j-^* t = (^1,^2, •••)• this moment, we temporally 
ignore the parameter z. Then the original KP theory states that a residue 
identity, called bilinear identity holds for (adjoint) wave functions: 

ResA w(z,t. A) u;*(z,t'. A) = 0, (4) 

where t' = {t'^.t^, ...) and the residue of A can be simply considered as the co- 
efficient of A~^ in the Laurent expansion. Since we have ignored the parameter 
z m. w and w* , the proof of (jl]) is the same as the original one given in [6]. 
Notice that the parameter z takes the same value for w and w* in (jlj). The 
key step in the proof of (jl]) is an important lemma. We recall it here [7]: 

Lemma 1. 

Resa P-Q* = Rcsa P(e«(*'^)) • Q(e"«(*'^)), (5) 

where P and Q are pseudo- differential operators. The residue with respect to 
d on the left hand side of ^ is defined as the coefficient of for a pseudo- 
differential operator. 

Actually, it is not difficult to prove that the bilinear identity (jl]) is equiva- 
lent to the KP hierarchy ( !Ta|) (see [7] for the proof). Namely, if we start with 
( lTa|) or ([2]), we can prove that wave functions defined by ([3]) will satisfy the 
bilinear identity dl]). Conversely, if we have wave and adjoint wave functions 
satisfying bilinear identity (|1]), then we can immediately know that the dress- 
ing operator W corresponding to the wave function w will satisfy (|2]), which 
makes the operator L defined by dressing form WdW~^ as a Lax operator for 
KP hierarchy. 

Then, if we consider KP hierarchy with squared eigenfunction symmetry, 
we can see that the dressing operator should satisfy another equation: 

d,W = qd-^rW. (6) 

In this case, we have: 
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Theorem 1. The KP hierarchy with a squared eigenf unction symmetry (J\) is 
equivalent to the following residue identities 

Resxw{z,t,X) ■ w*{z,t' , X) = 0, (7a) 

Resx Wz{z,t, \) ■ w*{z,t' , \) = q{z,t)r{z,t'), (7b) 

Resxw{z,t,X) ■ {q{z,t')w*{z,t',X)) = -q{z,t), (7c) 

Resxd-^ {r{z,t)w{z,t,X)) ■ w*{z,t' , X) = r{z,t'), (7d) 

where the inverse of d is understood as pseudo- differential operator acting on 
an exponential function, e.g., (rw) = (d~^rW)(e^^^'^^) . 

Proof. Frist, we will prove the residue identities (I7j) from ([T]), ([2]) and ([6]). 
Equation (17al) is easy to prove, since L satisfies the evolution equations of 
original KP hierarchy. 

To prove (]7b[) . we only need to show 

Res, w,{z, t. A) ■ (9,7 ■ • ■ dr,'w*{z, t. A)) = q{z, t) ■ • ■ d^^^riz, t) 

for arbitrary k > 1 and mj > 0. Here we have a key observation: the action of 
mixed partial derivatives 9^ ■ ■ ■ d^'' on w* can be written as Pmi---mkW*{z, t. A), 
and 

9,7 •• -97^^,1) = P„,...„, (r(^,t)), 

where Pmi---m^ is a differential operator in d. Then by noticing Wz{z,t,X) = 
qd~^rw{z,t, A) and Lemma [H we have 

ResA w;.(z,t,A)97---97w;*(z,t,A) 
=ResA q{z, t)9^V(z, t)w{z, t. A) • P„,...„^ty*(2;, t. A) 
=Resa g(2;,t)9"V(2,t)P^^...„^ 
=q{z,t)Pm,...mk{r{z,t)). 

Analogously, we have another residue identity: 

ResA w(z, t. A) wl{z, t'. A) = q{z, t)r{z, t'). (8) 

However this equation is equivalent to (iTbl) . so this equation is not included 
in (CD. 
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The residue identities (!7c|) and (17d|) can be easily derived from (17bp (and 
(IH])) if one substitutes Wz = qd~^rw and w;* = —rd~^qw* into ( iTbl) (or ([H])) 
and eliminate g or r. 

The proof for the reverse part in this theorem is written as the following 
proposition. □ 

Proposition 1. If there are functions q{z,t), r{z,t), 

w{z, t, A) = (1+$^ w,X-'j e«(*'^) and w*{z, t, A) = U+Yl J 

^ i>l ^ i>l ^ 

satisfying the identities ([^, t/ien t/iere exists a pseudo-differential operator W , 
and Lax operator L := WdW~^ such that L, q and r satisfy (Qp. 

Proof. Supposing W = 1 + X]i>i ""^j^"'' W = 1 + ^i>iW*d~\ it is easy to 
see w{z,t,X) = H/(e«(*'^)) and w*{z,t,X) = W {e~^^^^^'^) . Then from (JTa]) and 
Lemma [H we know the 9-residue Resg WW*d"^ vanishes for m > 1, which 
implies the negative part (W^W^*)_ is 0. Furthermore, the non-negative part 
of WW* is 1. So we have proved that WW* = 1, which means W = {W*)~^. 

Now, let's define L := WdW~^. We will prove Wt^ = —L'^W, which 
implies (fTa|) . From the definition form of W, it is easy to see that (Wt„ + 
L"iW)+ = 0. Again from ( 17al) we have 

ResaiWt^W-^ + L!^)„9'" = Resd{Wt^W-^ + {Wd''W-^)^)^d"' 
= ResxWt^e^ ■ {-d)'^{W*)-^e-^ + ResoiWd^'W-^ - {Wd''W-^)+)d'^ 
= ResxWt^e^ ■ i-dnw*)-'e-^ + ResxWd^'e^ ■ i-d)"'W*-^e-^ 
= Resxwt^{z,t, \){-d)'"w*{z,t, \) = 0, (for m > 0) 

so we obtain Wt„ = -L'^W. 

Next, we need to prove Wz = qd~^rW, which implies ( llbl) . It is easy to 
see that (WzW~^)+ = 0, and by using ( l7b| we have 

ResaWzW-^d"^ = ResAlV.e^ • (-9)™iy*-^e-« 
= ResAW,(^, t. A) • {-drw*{z, t. A) = g(-9)'"r (for m > 0), 

which means WzW~^ = J2m=oli~'^)"^i''")'^~"^~^ = qd~^r, namely, Wz = 
qd'^rW. 
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As the final step, we will prove ( !Tc|) and (lldp . We have already proved 
Wt„ = —L"iW, which implies Wt„ = L'^{w). Then according to (ITcl) . we find 



3 Bilinear Identity for the Extended KP Hi- 
erarchy 

In [23], we defined a new extended KP hierarchy. The key idea is to combine 
a specific k-th order fiow of KP hierarchy with the squared eigenfunction sym- 
metry fiow like f llbp . In this way, we found two types of KP hierarchy with 
self-consistent sources and their Lax representations. 

In [131 EH SO], the Hirota's bilinear equations for some (2+l)D integrable 
equations with sources are constructed by source generation procedure. For 
the case of KP equation, they found that these Hirota's bilinear equations 
correspond to the first and second type of KP equation with self-consistent 
sources. 

It is well-known that the bilinear identities are generating equations for 
Hirota bilinear equations of KP hierarchy. The natural questions are how to 
find the bilinear identities for the extended KP hierarchy and how to derive 
Hirota's bilinear equations from them. So in this section, we will give a de- 
tailed construction on how to derive the bilinear identities for the extended 
KP hierarchy (|9]). 

We remind that the extended KP hierarchy (for a fixed k) 



-ResA wt„(^,t. A) -d ^ {q{z,t')w*{z,t',X)) 
LI {-Resxw{z, t. A) ■ d^^ {q{z, t')w*{z, t'. A))) 
Ll{q{z,t)). 



In a similar way, we can find (fTd]) by using ( l7dl) . 



□ 



[Ll,L] {n^k) 
[LX + V, L] 
Ll{q) 



(9a) 
(9b) 
(9c) 
(9d) 
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is constructed by replacing an arbitrary fixed k-th flow dt,^ by dj^ where the in- 
flow is a combination of dt^.- and 9^-flow given by ( I9bl) . Functions q and r are 
eigenfunction and adjoint eigenfunction satisfying fl9c|) and fl9dl) . respectively. 
Note that q and r depend on rather than in the extendeded KP hierarchy 

©• 

Remark 1. /n order to simplify the notions, we will still use the symbols 
w{z,t,X), w*{z,t,X), q{z,t), r{z,t), L, W , etc., in this and the following 
sections, but they correspond to the case of the extended KP hierarchy for 
example, from now on, t = (^1,^2, • ■ • , tk-i,ik, tk+i, ■ • ■ ). 

To flnd the bilinear identity for the extended KP hierarchy ([9]), the key 
idea is supposing that L depends on the auxiliary variable z, whose evolution 
is given by ( ITb|) . 

Then we assume that the dressing operator for L with auxiliary variable z 
is given by = 1 + Xli>i t)9~'. And the evolution of W with respect 
to z and ik are given by (E]) and 

dt^W = -LtW + qd^^rW. (10) 

The definition forms for wave function and adjoint wave function are the same 
as ([3]) except for the functions depending on ik rather than tk e.g., ^(t,A) = 

Then we can prove the following theorem. 

Theorem 2. The bilinear identity for the extended KP hierarchy ^ is given 
by the following sets of residue identities with auxiliary variable z: 

Res\ w{z — ik,t, \) ■ w*{z — i'i^,t' , \) = 0, (11a) 
ResxW^{z-ik,t,X) ■ w*{z - i'^,t' , \) = q{z - ik,t)r{z - i'f.,t'), (lib) 
ResA w{z - tk, t, A) • {q{z - t^, t')w*{z - t'^, t', A)) = -q{z - 4, t), (11c) 
ResA d-' {r{z - h, t)w{z - h, t, A)) ■ w*{z - i'k, t', A) = r{z - i',, t'), (lid) 

where the inverse of d is understood as pseudo-differential operator acting on 
an exponential function, e.g., d~^{rw) = {d~^rW){e^) . 
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e-«(t,A) 



Proof. We notice that -^w{z — ik, t, A) = Wtf,{z — ik, t, A) — Wz{z — t^, t, A) = 
L'!l{w). The A-residue of [-^w{z — ik,t, X)]w*{z — tk,t,X) simply vanishes 
according to Lemma [TJ The same is true for the A-residues of arbitrary mixed 
derivatives [ ^^mi^jila w{z — tk, t, X)]w*{z — ik, t, A). Then identity (11 lap holds. 

The proofs of flllbl) - fllldl) are almost the same as the proofs of fl7b|) - fl7d|) . 
One should notice again that — tk, t, A) = □ 

Theorem 3. Suppose that there are (adjoint) wave functions 
w{z,t,X) = (Y,w,{z,t)xAe^^'''\ w*{z,t,X) = (Y,w,{z,t)xA 

^ i>0 ^ ^ i>0 ^ 

(with Wq = Wq = 1), q{z,t) and r{z,t) satisfying the bilinear relations / fli]) . 
then the pseudo-differential operator L = WdW~^, W = '^i>o'Wi{z,t)d~^ , 
and functions q and r give a solution to the extended KP hierarchy 

Proof. Let us define W := ^iyQWi{z,t)d~'^ . With the same argument as in 
Theorem [U, we will find that W = {W"^)*. Then from (lllbl) and Lemma [H 
we find that for any m > 0, 

ResaW.W-^-d)'^ = g9'"(r), 

which means (W.W-^)^ = qd~^r. Notice that the non-negative part of 
WzW~'^ vanishes, so we have Wz = qd~^rW. 

Next, from the coefficient of Taylor expansion of (11 lap , we find 

ResA -^w{z - ik, t. A) • d'^w*{z - ik, t. A) = 0. 
dtk 

d - d - 

By reahzing that ---^w{z-tk,t,X) is {-^W {z -tk,t) + L^W) ex^i^i^it, X)) and 

dtk dtk 
using Lemma [U we find 

Refia{^W{z - ik, t) + LtW)W-\-d)'^ = 0, (for m > 0) 
dtk 

which means -^W{z-ik,t) = -L^_W. Hence dt^W{z,i) = {-L'l + qd~^r)W . 

dtk ~— . 
For equations (jHc]) and (l9d|) . the proof can be done by differentiating (IllcD 

and (llldp directly. Since dt^w = L'^iw), dt„w* = -(L")*(w*), the rest of 

proof is obvious. □ 
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4 Tau-Function for the Extended KP Hierar- 
chy 

The existence of r-function for original bilinear identity of KP hierarchy is 
proved in [6]. In our case, the wave function w{z — ik, t, A) and w*{z — ik, t, A) 
satisfy exactly the same bilinear identity ( II lap as the original KP case if one 
considers z as an additional parameter. So it is reasonable to assume the 
existence of r-function and make the following ansatz: 

Hz - 4, t, A) = ^^'~'; + ^'V ^^^^ • (12a) 

t{z - tk,t) 

t(z - h - TTk^t + \X\) 

w*{z - t,, t, A) = ^ ' f ' , ^ ■ exp(-e(t, A)), (12b) 

r{z - tk,t) 

where [A] = iP"' sP"' ' ' ' ) ■ According to we should make further as- 
sumptions: 

q{z, t = r{z, t = 12c 

Then, similar with [1], we have the following results. 

d-\r{z - 4, tMz - t„ t. A)) = ^iiZii±i4^e^(M)^ (12d) 

\t{z - tk,t) 

d-\q{z - 4, t)w*iz - 4, t. A)) = Z^^l^lpm^^Me-^i^^). (I2e) 

At{z - tk,t) 

To find Hirota's bilinear equations for extended KP hierarchy (Q, we sub- 
stitute ([I2D into ([n]), and get 

ResA f{z,t- [A]) f{z, t' + [A]) e«(t-t''^) = 0, (13a) 
ResA f,{z,t- [A]) f{z,t' + [A]) e^^*"*''^) 

-ResAf(z,t-[A]) (a,logf(^,t)) f{z,t'+[X]) e^^'^''''^ 

= a{z,t) p{z,t'), (13b) 

ResA X~^f{z, t - [A]) a{z, t' + [A]) e«(*-*''^) = a{z, t) f{z, t'), (13c) 

ResA X-'p{z,t- [A]) f{z,t' + [A]) e«(*"*'^) = p{z,t') f{z,t). (13d) 
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Here the bar over a function f{z,t) is defined as f{z,t) = f{z — tfc,t), e.g, 
f(z,t-[A]) =r(z-(4-^),t-[A]),f(^,t' + [A]) =r(^-(4 + ^),t' + [A]). 

After setting t as t + y and t' as t — y in f fT3|) with y = {yi,y2, ■ ■ ■), We 
can write ( !T3|) as the following systems with Hirota bilinear derivatives D and 
A's: 

5^P.(2y)p.+i(-^) exp [ J2 VjDj ) f{z, t) • fiz, t) = 0, (14a) 
i>o \i>i / 

^Pi(2y)pi+i(-£))exp i^yjDj J f^(2;,t) ■ f(2;,t) 
i>0 \j>i J 

- 5^p.(2y) {d, log f (2;, t + y))p,+i{-D)f(z, t + y) • f (z, t - y) 

j>0 

= exp yjD^ a{z, t) • p(2, t) (14b) 
^Pi{2y)pi{-D)exp I ^yjDj j f(2;,t) • o-(2;,t) = exp I ^yjDj j o-(2;,t) ■ f{z,t) 

i>0 \j>l J \j>l J 

(14c) 

^Pi{2y)pi{-D)exp I ^yjDj J p(2;,t) ■r(2;,t) = exp I ^yjDj J f(2,t) ■p(2;,t) 
i>0 \i>i / \i>i / 

(14d) 

where D = (Di, • ■ • ), -Dj is the well-known Hirota bilinear derivative 

Dif ■ g = ftiQ — fOu, and Pi{y) is the z-th Schur polynomial, whose generating 
function is given by 

00 00 
exp^y.A^ = ^Pi(y)A\ 

i=l i=0 

Remark 2. The zero-th order term in {Ijb^ (with yj = y j) can be written 
in the following forms with Hirota 's operator 

1 

ap + D^T^ .T = ap + D^T^, ■ r = ap + -D.^D^t ■ r = 0. 



Note that (1141) is the generating equations of Hirota's bilinear equations for 
the extended KP hierarchy (jO]), where the dependence of auxiliary parameter 
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z can be eliminated if one wants to convert bilinear equations to PDEs. In the 
following sections we will give several examples to show the Hirota's bilinear 
forms for the nonlinear evolution equations in extended KP hierarchy and how 
it can be transformed back to nonlinear PDEs (which will correspond to the 
well-known KPSCS-I and 11, etc.). 



Example 1 (First type of KP equation with a source (KPSCS-I) [T] 
[32] . i.e., the extended KP hierarchy ([9]) for n = 2 and /c = 3). The Hirota 
equations for the KPSCS-I can be obtained as 

DxTz ■ T + crp = 0, from (TJ^ with yj = 

(D^ + 3Dl - ADx{D-t, - D,))t • r = 0, from ^ m 

{Dt^ + DI)t • ct = 0, from \14c\ ) in y2 

{Dt^ + Dl)p ■ r = 0, from in y2 

where is Hirota's derivative, i.e., Dzf{z) ■ g{z) = fzQ — fdz- Note that from 
the definition of f, we know that Dp^f ■ f = {Di^ — Dz)t ■ t, which interprets 
the appearance of this term in the second equation. 



Example 2 (Second type of KP with a source (KPSCS-II) [E 
i.e., the extended KP hierarchy ([9]) for n = 3 and k = 2). The Hirota equations 
for the KPSCS-II can be obtained as 

DxTz ■ T + ap = 0, from Ijl^b]} with yj = 

(D^ + 3(A-, - DzY - AD.,Dt,)T • r = 0, from ([1^) m y^ 

{{Dt, - Dz) + DI)t -(7 = 0, from m y2 

{{Dt, - Dz) + Dl)p • r = 0. from (IJ^ m y2 

(4^3 -Dl + 3D,.(A-, - Dz))t ■a = 0, from (TJ^ in y^ 

{ADt, -Dl + 3DxiDt, - Dz))p • r = 0. from m 

It seems that the Hirota bilinear equations obtained here for KPSCS-II are in 
a simpler form comparing with the results by Hu and Wang fl^ . 
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Example 3 (The extended KP hierarchy (Q ioi n = 4 and k 
Hirota bilinear equations for this system can be obtained as 



2 [25]). The 



DxTz • r + crp = 0, 



T ■ T 



0, 



3D,Dt, - 2{Di, - D,)Dt, - Dl{Dt, - 

4^3 -Dl + 3D,(A-. - D,)) r ■ a = 0, 
4A3 -Dl + 3D,{Dt, - D,)) p • r = 0, 
I8A4 + Dt- QDliD-t, - Dz) + 3(A-. - 



T ■ T 



SD^Dt,) p-T - 



0, 
0. 



from |j^6[ ) with yj = 
from (IJa) in y^ 
from (14^ in y^ 
from \i4c\j in y2 
from (14^ in y2 
from il4c\) in y^ 
from { 14d\ ) in y^ 
from ( [i^c[ j in 
from imd^ ) in 



5 Back to nonlinear equations from Hirota's 
bilinear equations 

To see the nonhnear PDEs corresponding to the Hirota's bihnear equations in 
the previous examples in last section, we convert the bilinear equations back 
to nonlinear PDEs in this section. It is not as simple as generating special 
bilinar equations. We follow the way given in [12]. 
Consider the following identities 

exp j P-T = e2™«h(E.^«^0 . e^«^'^> (p/r) , (15a) 

cosh (^6,D}j T-T = e2=°«h(E.'5.aoiogr_ ^^51^^ 



dcf 



Remark 3. Note that in this section, Dif ■ g = ft^g — fgi^- 

The identities (fT5l) are easy to prove. One should notice the definitions of 
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operator Di, for example, f llSaP is proved by checking from left hand side: 
l.h.s. = p{ti + 5i,t2 + 62, - ■ ■ - 5i,t2 - 62,- ■ 

p(ti + 61, t2 + (^2, • ■ • ) 



r.h.s. = exp [(e^'''^^' + e ^t^^^') logr] 

= exp [logr(ti + 5i,t2 + ^2, ■ • ■) H 
p(ti + 61, t2 + ^2, ■ • ■) 



T{tl+5i,t2 + 52,---) 

logr(ti - 5i,t2 - ^2, • ■ •)] 



Titi + Si,t2 + S2,---) 

= p{ti + 5i,t2 + 62, - ■ ■ )r{ti - 5i,t2 - 62,- ■ ■)■ 
Using the transformations u '= 9^ logr {x = ti), r =^ p/r, q '= cr/r, we 



can rewrite fllSp into the following form. 
-2!^ p-r = exp 



n=0 

00 



00 v-^ / r- r% \ 2n 



2n! 



n=0 



2n\ 



\2n 



-T ■ T = exp 



n=0 



n=0 



2n\ 



(16a) 
(16b) 



Thus, relations between the Hirota bilinear derivatives and usual partial deriva- 
tives with respect to original variables are established. We list some of these 
relations as follows. 



D'fT-T 




D.D^r ■ r 


DiD^T ■ T 




r2 


r2 


r2 


2u 


2d~^U2 




29-1^4 


D2D^T ■ T 


DiD2T ■ T 


D^r -r 






t2 


t2 






2tii,2 + 12M(a-^M2) 


2mi,i + 12^2 




D,p-T 


i>'2P ■ r 




D^p ■ r 






t2 


r2 


ri 


^^2 




r4 


Dip-T 


/^i/^2P ■ r 




U^p ■ T 




r2 


t2 


r2 


Til + 2ur 




ri^3 + 2r5"^M3 


r2,2 + 2rd~'^U2,2 


D'tp-r 


DID2P ■ T 








t2 


t2 






?'i,i,2 + 2ur2 + Arid^^U2 


ri,i,i,i + 12Mri,i 








2(6^2 + Mi,i)r 
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where the subscript i, j ■ ■ ■ of Mjj-,... "^id,--- denote the derivative with respect 
to the corresponding variables tj, tj • ■ • . 

Example 4. ExampleUl can be translated back to nonlinear PDEs as (y = t2) 

dzd~^u + gr = 0, 

{uxxx + 12uux' - 4uf-3 + Auz).j, + "iuyy = 0, 
Qy = Qxx + 2Mg, 

Ty "l^xx ^UT. 

If we eliminate auxiliary variable z from the above equations, we get the first 
type of KP equation with a source m\ [13, [23. [2^ 

(4^43 - u^^^ - 12nn^)a; - 'iuyy + 4:{qr)^^ = 0, 
qy = qxx + 2Mg, 

Ty T^xx '2,UT^ 

whose soliton solutions can be obtained by dressing method JMjj and Hirota 
method J77| /- 

Example 5. ExamplelE can be translated back to nonlinear PDEs as (y = 

t2,t = t3) 

dzd~^u + = 0, 

{uxxx + I2uux - ^ut^)x + 'i{dy - dzfu = 0, 
Qy-Qz = Qxx + 2Mg, 

' y ' z ' XX ^u,/ , 

- Aqt + g^xx + Quq.^ + Sg^y - Sg^,^ + Qqdyd'^u - Gqd^d'^u = 0, 
— r^^x — 6ur^ + 3r^y — Sr^^ + Grdyd^^u — Grd^d^^u = 0. 

By eliminating auxiliary variable z, we reach the second type of KP equation 
with a source (T^IMIMIM 

{4:Ut - u^xx - I2uu^)x - Suyy = 3 [(gr)y + (g^^r - gr^^Ola; ' 
qt = Qxxx + ^uq^ + '^u^q + ^qd'^Uy + ^gV, 
n = r^xx + 3ur^ + ^u^r - ^rd~^Uy - ^r^g, 
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whose soliton solutions can he obtained by dressing method 12^ and Hirota 
method JT^ . 

Example 6. Example\E can be translated to nonlinear PDEs as (y = t2) 

dzd~^u + qr = 0, 

2u^^ + 12^2 + 6d-'^[{dy - d,fu] - 8d-^Ut., = 0, 

6d"^ut^ - 49~^(Mt3,y - ^t3,z) - I2ud^^{uy - u^) - 2{u^y - u^z) = 0, 

qy = qz + qxx + 2ng, 

y ' z 'XX ^ u,/ , 

4gt3 = qxxx + Quq^ + 3 [q^y - q^^^ + 2qd~^{uy - u^)) , 
"irt^ = r^xx + Gur^ - 3 {r^y - r^^ + 2rd~^{uy - u^)) , 
18gt4 = Qxxxx + 12Mgxx + 2Aq^d~^{uy - u^) + (12m^ + 2u^^ + IQd'^ut:, 

+ 6d~'^{dy - dzYu)q + 12u{qy - q^) + Q{qxxy - Qxxz) + ^{dy - d^fq + Sg^tg, 
18rf4 = -r^xzx - 12ur2;2; + 24r^(9"^(uy - u^) - [l2u^ + 2n^^ + IQd^^'at,, 

+ 69"2(9^ - d,Yu)r + 12M(ry - r,) + 6(r,.,.^ - r,.,.,) - ?>{dy - d^fr - 8r,t, 

Eliminating the dz and dt^ terms, we get the extended KP hierarchy ^ for 
n = 4 and k = 2 IE 



2Ut^ =d ^Uyyy + U^^y + 4(M^)j^ + 4(9 ^ Uy')Ux + 8{uqr)x + d ^{qf')yy 

+ 2(ga;a;r + qrxx)x + {rq^ - qrx)y, (17a) 
qu =Qxxxx + "^uq^x + (4ux + 4gr + 2d~^Uy)qx + {d~^Uyy + d~^{qr)y 

+ Au^ + 2m^.^. + Uy)q, (17b) 
-^^44 ='^xxxx + 4Mr^x + (4mx - 4gr - 2d^^Uy)rx + {d~'^Uyy + (9"^(gr)j^ 

+ 4^2 + 2n^.^. - nj,)r, (17c) 

whose Wronskian type solution (including soliton solutions) can be obtained 
by dressing method I2^ . The correctness of can be verified from the zero- 
curvature equation of extended KP hierarchy ^ with n = A and k = 2: 

B^,i, - + [^4, + qd-\]+ = 0, (18a) 
duq = B^{q), (18b) 
dt,r = -Bl{r). (18c) 



17 



In the expressions of [T^) . if we eliminate U2 and from term by using 
and d terms, we can verify that fT8\} is exactly {11}. 



6 Conclusion and Discussions 

In this paper, we constructed the bihnear identities f|TT]) for the extended KP 
hierarchy ([9]) defined in [25] . The bihnear identities f lTTj) are used to derive the 
Hirota's bihnear equations f|T^ for ah the zero- curvature forms of the extended 
KP hierarchy ([9]). We have shown that the Hirota's bihnear forms in Example 
[U and Example [2] correspond to the KP equation with a self-consistent source 
(KPSCS-I and II). After translating the Hirota equations back to nonlinear 
PDEs, the correctness of our bilinear forms (fTT!) are verified. Another forms of 
Hirota's bilinear equations for KPSCS-I and KPSCS-II have been given in [T3] 
by using source generation procedure., where the auxiliary functions Pi and Qi 
are introduced. In this paper, by introducing an auxiliary fiow (S^-fiow), the 
bilinear identity for the whole extended KP hierarchy and a simpler Hirota's 
form are obtained. To show the validity of our method, we gave an extra 
example in the extended KP hierarchy (Example [3j). 

There are some important applications for bilinear identities of extended 
KP hierarchy. As we know, the quasi-periodic solutions for KP hierarchy can 
be constructed by using method in algebraic geometry, where the construction 
of wave functions (or Baker- Akhiezer functions as in quasi-periodic cases) are 
intimately related with bilinear identities, Riemann surfaces and divisors on it. 
It is very interesting to consider the quasi-periodic solutions for the extended 
KP hierarchy ([9]) when bilinear identities have been obtained in this paper. 
Another interesting problem is to consider the bilinear identities for other 
extended hierarchies, such as BKP, CKP, 2D Toda and discrete KP, etc. We 
would like to investigate these problems in the future. 
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